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ABSTRACT 

The one-dimensional contact process is analyzed by a cluster approximation. In this 
approach, the hierarchy of rate equations for the densities of finite length empty in- 
tervals are truncated under the assumption that adjacent intervals are not correlated. 
This assumption yields a first order phase transition from an active state to the ad- 
sorbing state. Despite the apparent failure of this approximation in describing the 
critical behavior, our approach provides an accurate description of the steady state 
properties for a significant range of desorption rates. Moreover, the resulting criti- 
cal exponents are closer to the simulation values in comparison with site mean-field 
theory. 



The contact process (CP) is an irreversible lattice model involving nearest neigh- 
bor interactions only This model incorporates spontaneous desorption and 
nearest-neighbor induced adsorption. This stochastic process can be used to mimic 
epidemic spread as well as catalytic reactions. This model belongs to a general class 
of nonequilibrium models exhibiting a continuous phase transition. Near the criti- 
cal point, the system exhibits divergence of spatial and temporal correlations. Such 
properties, conveniently characterized by critical exponents, can be used to classify 
different models. The CP belongs to the same universality class as Sclogel's first 
model [[|, Reggeon field theory ||, directed percolation ||, and the ZGB model || 
of catalysis. Field theoretic renormalization group studies J7], |8|] provide considerable 
understanding of the critical behavior of the CP. However, the best estimates for the 
characteristic exponents were found numerically by Monte Carlo simulations M and 
by series expansion analysis [[H], [Tl] • 

Motivated by the incomplete theoretical understanding, we introduce an approx- 
imate approach to the CP. We study the temporal evolution of the density of empty 
intervals. The corresponding rate equations lead to an infinite hierarchy of equations. 
By writing the density of pairs of neighboring empty intervals as a product over single 
interval densities, we obtain a closed set of equations. We use the generating function 
technique to obtain the steady-state properties of the system. Within this approxi- 
mation, the system exhibits a discontinuous phase transition from an active state to 
the empty state. As the system approaches the critical point, the relaxation time, 
associated with the temporal approach to the final state, diverges. Consequently, at 
the critical point, an anomalously slow decay towards the final state takes place. We 
find the corresponding kinetic exponent by scaling techniques, as well as by numeric 
integration of the rate equations. 

We compare the cluster approximation predictions with the results of site mean- 
field theory and with series analysis of this process. Despite the failure to predict a 
continuous transition, the cluster approximation provides a good approximation for 
the final density and the empty interval density for a reasonable range of desorption 
rates. Moreover, the resulting estimates for the critical exponents are closer to the 
numeric values in comparison with site mean-field theory. The cluster approach is 
also applicable to generalizations of the contact process, such as the A model and 
the N3 model. We verify that the resulting critical behavior of these processes is 
identical with the contact process. Our approach is advantageous since it can be 
improved systematically by considering the evolution of higher order empty interval 



1 



• ° • • • o 

Figure 1: The Contact Process. 

densities. 

In the CP, a particle desorbs spontaneously with rate A. On the other hand, a 
particle adsorbs at a given site at a rate proportional to the number of neighboring 
occupied sites. In other words, the adsorption rate at a particular site is given by 
n p /n s , with n s the total number of neighboring sites and n p the number of neigh- 
boring particles. Since every neighboring particle contributes independently to the 
adsorption rate, this stochastic process can viewed as an interacting particle system 
with nearest-neighbor interactions only. The above process possesses an adsorbing 
"vacuum" state: once the system reaches the empty state, adsorption becomes im- 
possible. In sufficiently high dimensions, neighboring sites are not correlated, and 
the density follows from dp/dt = p(l — p) — Xp. The adsorption term represents the 
density of vacant sites that neighbor an occupied site. This site mean-field theory 
(SMF) gives a steady state concentration equal to 1 — A. Hence, at A c = 1 this 
process undergoes a simple continuous transition. 

We consider the one-dimensional case only (see Figure [I]), where studying the 
density of empty interval has proven useful in adsorption processes [Q, as well as 
in reaction processes JH|]. Denote by E n (t) the probability that a randomly chosen 
string of n sites is empty (see Figure H). We emphasize the fact that the actual empty 
string might be of length larger than n. Let us also consider R n {t), the probability 
that a random string of length n+1 has n consecutive vacant sites and a particle at 
the extreme right site. For symmetric processes, such as the CP, R n (t) also represent 
the probability of finding an empty string with a particle at the extreme left. These 
two interval densities are related by 

R n = E n — E n+ i or E n = 1 — Rq Rn-x, (1) 

for n > 0. For n = 0, the definition of E n is trivial, leading to the following conditions 
satisfied by the empty interval densities, 

oo 

E = l or J2 R n = 1 - ( 2 ) 

n=0 
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E n = Prob( Q • • • Q ) R n = Prob( Q ■ ■ • Q # ) 

n n 
Figure 2: Empty Interval Densities 

The condition for R n is obtained by using Eq. ([!]) and noting that the sum over E n 
reduces to an alternating series. The above interval densities are useful in describing 
macroscopic properties; for example, the concentration is given by p(t) = Ro(t). 

We write the rate equations governing E n (t) by considering the adsorption and 
the desorption processes separately. Adsorption can contribute only to loss of empty 
intervals. Empty intervals of length n can be destroyed when a particle is desorbed at 
the edge of the interval. This occured only where the empty interval has an occupied 
site at its edge. By taking into account contributions from adsorption at either 
boundary, we find cLE n /<it| ads = —R n . Desorption, on the other hand, leads only 
to creation of empty intervals. When a particle desorbs, the two empty intervals it 
borders create a larger empty interval. We define the pair density E^ m as the density 
of two neighboring empty intervals, with lengths at least equal to / and m, separated 
by a single particle. Note that from this definition the relation E 0tU = E nt0 = R n is 
satisfied. In terms of this pair density, the increase in the density of empty intervals of 
length at least n due to desorption is described by dE n /dt\ des = \ Yd=o -S^n-i-i- We 
then approximate the pair density, E^ m , by the product E^ m = RiR m / Ro, where the 
factor 1/Ro ensures the normalization condition Eq^ = R n . Combining contributions 
from both adsorption and desorption yields the following rate equation for the empty 
interval density, 

dE A n ^ 

-£ = -Rn + -B-Y. RiRn-i-h n > 0. (3) 
at Kq l=0 

The steady-state properties can be obtained by requiring that the time derivative 
in Eq. @ vanishes. We introduce the generating function R(z) = J2 n RnZ n /Ro, with 
Rn being the steady-state interval densities. Dividing Eq. (|3]) by R , summing over 
all n, and solving the resulting quadratic equation yields 



AM = '"^f^ . » < V4. (4) 

The normalization condition of Eq. (|2]) shows that p = Rq = 1/ R(z)\ z=1 and, conse- 
quently, the concentration is given by 

l + v / T^4A ^ 1/A 
P = o > A < V 4 - (5) 
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Figure 3: The cluster approximation density (solid) vs. the actual CP density 
(dashed). The latter density represent the [12,12] Pade approximant obtained from 
the perturbation study of ref fil l. 



For A > A c = 1/4 the system indeed exhibits a transition to the absorbing state, 
where the concentration vanishes. The nature of this transition is discontinuous, 
while for the actual CP, the transition is continuous. For the CP, the approach to 
the critical density is an algebraic one, 

p- Pc ~(\ c -\f. (6) 

Note that while for a continuous transition to a vacuum state one has p c = 0, the 
critical density of the cluster approximation is finite, p c = 1/2. Extensive power 



series studies suggest (3 = 0.277, and A c = 0.3032 JR], [TT|. The corresponding values 
obtained by the cluster approximation, (3 = 1/2 and A c = 1/4, are closer than the 
SMF values (3 = X c = 1. While the critical point depends heavily on the microscopic 
definition of the process, the critical exponents are universal. Applying the cluster 
approximation to variants of the CP always yields (3 = 1/2. 

In Figure [| we plot the cluster approximation density vs. the series study density. 
For desorption rates < 0.2 both curves are practically identical and, for example, at 
A = 0.2 the relative difference is less than 0.5%. We conclude that despite the 
failure near the transition point, the cluster approximation is useful in describing the 
process for a substantial dynamic range. Another way to determine the accuracy of 
the approximation is by expanding Eq. ([5|) as a power series in A n , and comparing 
to the coefficients, obtained by the series expansion study 0. We find that both 
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expansions are identical to the third order, 

J 1 - A - A 2 - 2A 3 - 5A 4 + 0(A 5 ) CA, 
p(A) = \l-A-A 2 -2A 3 -4|A 4 + 0(A 5 ) CP. " ' 

Since R n is of order A n , we expect similar correspondence between the approximate 
interval density from our rate equations and the actual contact process densities. To 
test this, we have performed a Monte-Carlo simulation with 10 5 particles at A = 0.2. 
Indeed, the average over 200 different realizations yields values for R n that agree 
with the approximate density to within 0.5%, for n = 0,1, 2, 3. The quality of the 
approximation gradually decreases as n increases and for example for n = 4 the 
discrepancy is 4%. 

The interval density also exhibits an interesting critical behavior. By expanding 
the generating function of Eq. one finds R n = pX n (2n)\/n\(n + 1)! for A < A c . 
This expression shows that for A < A c the interval density depends exponentially on 
the interval length, 

R n ~ n^ /2 exp(-4An) A < 1, (8) 

where A = A c — A. The Stirling formula was used to obtain the above form of the 
interval density. At the critical point, power-law decay is recovered, R n ~ n~ 3 ^ 2 . 

These static properties are closely related to the kinetics of the system. Far from 
the critical point, the density rapidly relaxes to its steady-state value from any initial 
conditions. As the system approaches the critical point, the relaxation time diverges 
and a power-law decay of the concentration takes place. We thus expect that for 
sufficiently large times and close enough to the critical point, the temporal approach 



to the final state is given by the scaling form pO, 14 



p-p^r'MAt 1 '")- (9) 

In other words, the critical exponent 5 characterizes the critical kinetics, while the 
exponent v characterizes the sub-critical relaxation time. The above scaling form 
should match the steady-state form of Eq. @ at large times and hence we conclude 
that ip(x) ~ x 13 for x 3> 1. To cancel the temporal dependence, the scaling relation 
j3 = 5u must be satisfied. Furthermore, Eq. (H) indicates that the relaxation length 
associated with the steady-state interval density diverges as 1/A as the system ap- 
proaches the critical point. Thus, it is natural to assume that R n (t) depends on time 
through a rescaled size, n — > nt~ a , as well as a rescaled adsorption rate. Noting that 
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Figure 4: The critical approach towards the steady-state, 
rate equations (dots) for the case A = A c = 1/4 is plotted, 
slope —1/3 is also shown (solid). 



Numeric solution to the 
For comparison a line of 



Ro = p, we postulate the following scaling behavior for R n (t): 

R n (t) ~ <j)(At 1/u ,nt- a ). (10) 



On the other hand, Eq. (|8|) indicates that the steady-state density depends on the 
size and the rate only through the variable nA. Hence, <fr(x,y) ~ <p(xy), and by 
eliminating the temporal dependence we find the scaling relation olv — 1. 

Thus far, our scaling analysis involved matching the anticipated kinetic behav- 
ior to the exact steady-state properties. To determine the critical temporal decay, 
P — Pc ~ t~ 5 , we study the rate equations at A c = 1/4. Using the duality relations 
between E n and R n (see Eq. (JJ)), the rate equations for A = A c can be rewritten in 
terms of R n only, 



n— 1 ^-p^ n— 1 

^ = TR~ ^ R l R n-l~l ~ R n- (11) 

To analyze this equation by scaling techniques, we match the leading asymptotic 
terms in both sides of the above equation. The left hand side is governed by n terms 
of order R n . Hence, by taking into account the time derivative, we conclude that 
the left hand side is proportional to nR n /t. The right hand side is dominated by 
the first few terms in the expansion, namely / n and n — I <C n. Therefore, we 
approximate the sum by 2i?„(_R + Ri + • • ')/ Rq — 2R n /Ro, using the normalization 
condition of Eq. (Q). Finally, we write the resulting expression R n (1 — l/2R (t)) in 
terms of the concentration, 

Ti 

-R n oc R n (p - p c ). (12) 
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V 


CP 


0.277 


0.160 


1.735 


CA 


1/2 


1/3 


3/2 


SMF 


1 


1 


1 



Table 1: Static and Critical decay exponents obtained by series studies (CP), Cluster 
Approximations (CA) and Site Mean-Field (SMF). 

We conclude that nt^ 1 ~ ~ t~ s , or equivalently a + 5 = 1. The three scaling 
relations yield the following exponents, 5 = 1/3, v = 3/2 and a = 2/3. Numerical 
integration of the rate equation at A = A c confirms the scaling prediction for 5 
(see Figure |]). In Table |I|, we compare the exponents, that result from our cluster 
approximation (CA), with the corresponding series expansion and the SMF values. 
We conclude that the CA exponents provide significant improvement in comparison 
with those from SMF. 

Recently, several variants of the CP where introduced and it was shown that 
they belong to the same universality class as the original model. In these models, 
the adsorption process is modified while the desorption process remains unchanged. 
The adsorption rate is set to 6/2 if only one neighboring site is occupied, while 
the rate remains unity in the case when both neighboring sites are occupied. The 
empty interval method can be easily generalized to this case and we merely quote 
the resulting subcritical steady state density 



P 



e + Je 2 - 4\(6 - t]) 



with 



V 



A(l-0)(0-2A) 



(13) 



2(6-t]) ' ' 6 2 + (I - 8)(6 - 2\y 

For the basic CP, 9 = 1, we recover Eq. (j5[). Clearly, the discontinuous nature of the 
transition is independent of the microscopic details of the model and the exponent 
(3 = 1/2 is indeed robust. The critical point, however, depends on 6 and can be 
found by equating the square root in Eq. ( |13"D to zero. If the adsorption rates are 
independent of the number of neighboring particles, 6 = 2, the predicted critical 
point is A c = 0.4608, while series studies yield A c = 0.574 for this so-called A model. 
In the case where 6 = 1/2 (the N3 model), the critical point is A c = 0.1366, while 
series studies yield A c = 0.162. For the general 6 case, the cluster approach yields 
less accurate estimates for the concentration than for basic CP. 

In summary, we have presented an approximate approach to the contact process. 
Based on the assumption that neighboring empty intervals are not correlated, we 
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solved for the steady state properties. In addition, the kinetic approach towards the 
steady state was found by scaling techniques. The above approximation is valid for 
a significant subcritical range. The cluster approximation predicts a discontinuous 
transition but gives improved exponents in comparison with simple site mean-field 
theory. 

The cluster approximation can be systematically improved by considering higher 
order interval densities. Indeed, our preliminary results indicate that the second 
order cluster approximation yields a continuous phase transition. Moreover, the 
approximation appears to be valid over a significantly larger range of desorption 
rates. 
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